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Abstract 

A magnetization equation for a system of spins evolving non-adiabatically 
and out of equilibrium is derived without specifying the internal interactions. 
For relaxation processes, this equation provides a general form of magnetization 
damping. A special case of the spin-spin exchange interaction is considered. 
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1 Introduction 



Magnetization dynamics is conventionally described by the Landau-Lifshitz-Gilbert 
(LLG) equation [1],[2], which provides a plausible phenomenological model for many 
experimental results. Recently, the LLG equation and the Gilbert damping term have 
been derived from an effective Hamiltonian including the radiation-spin interaction 
(RSI) [3]. It has been assumed there that the spin system maintains quasi-adiabatic 
evolution. 

Various kinds of relaxation processes are usually melded together into a single 
damping term. Relativistic relaxation processes result in the Gilbert damping term 
with one damping parameter, while for the case of both exchange and relativistic 
relaxation the damping term is a tensor with several damping parameters [4] , [5] . 

The relaxation processes are specified by interactions of spins with each other 
and with other constituents of the magnetic system. A derivation of the damping 
term from first principles should therefore start with a microscopic description of the 
interactions. Even though such microscopic derivation of damping has been performed 
for some relaxation processes (for instance, [6]), a full version of derivation for the 
Gilbert damping term has not yet been given, and in particular for a system in non- 
equilibrium. 

In the present talk, we aim to derive a magnetization equation for a general non- 
equilibrium spin system without specifying the interaction Hamiltonian and related 
relaxation processes. We start with a system of spins precessing in the effective 

— * 

magnetic field H e ff neglecting for a moment mutual interactions. Then at a fixed time 
later interactions in the system are switched on and influence the original precessional 
motion. 

The interactions are assumed to be time- dependent, and the spin system evolves 
non-adiabatically out of equilibrium trying to relax to a new equilibrium magneti- 
zation. We perform a transformation, which is analogous to the one used in the 
transition to the interaction picture, to connect the density and magnetic moment 
operators before and after the time, when a new non-equilibrium dynamics starts, and 
to find an explicit expression for the interaction contribution to the magnetization 
equation. 

2 Magnetization Equation 

Let us consider a quantum spin system defined by 

H = H + \H I , (1) 

where Ho is the Zeeman Hamiltonian describing the interaction of spins with an 
effective magnetic field 

Wo = -7E&'4//(*), ( 2 ) 
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7 being the gyromagnetic ratio and Si being the spin operator of the ith atom. The 
effective magnetic field H e ff is given by the energy variational with magnetization , 
Hg/j = —5E(M)/5~M., where E(M) is the free energy of the magnetic system. This 
field includes the exchange field, the anisotropy field, and the demagnetizing field, as 
well as the external field, H e:rt . 

The Hamiltonian Hi represents other possible types of interactions, which can in- 
clude, for instance, higher order spin-spin interactions. The interaction terms included 
in Hi are in general time-dependent, being switched on adiabatically or instantly at 
a fixed time t . The parameter A in (1) can be chosen small in order to take into 
account the higher order effects perturbatively. 

We introduce next the magnetic moment operator 

M ss (3) 

which is the response of the spin system to the external field. The magnetization is 
defined as an ensemble average of the response 

M=(M) = ^Tr{pM}, (4) 

where V is the volume of the system, while p is the density operator satisfying the 
quantum Liouville-von Neumann (LvN) equation 

ih^ + [p,H]_ = 0. (5) 

For systems in equilibrium, the Hamiltonian itself satisfies the LvN equation and 
the density operator is expressed in terms of the Hamiltonian. For non-equilibrium 
systems, the density operator is constructed by making use of the time-dependent 
adiabatic invariants [7], [8]. 

To derive the magnetization equation, we proceed as follows. We perform, on p(t), 
the transformation 

p^ Pmt = U(t ,t)p(t)U(t,t ) (6) 

defined by the operator 

U(t ,t)=Texp{^ f'drHoir)}, (7) 
n J t t) 

where T denotes the time-ordering operator. For systems with Ho constant in time, 
the operator U(to,t) = exp{(i/h)Ho(t — t )} leads to the interaction picture, which 
proves to be very useful for all forms of interactions since it distinguishes among the 
interaction times. For our system with both Ho and Hi dependent on time, the 
operator (7) plays the same role, removing the unperturbed part of the Hamiltonian 
from the LvN equation. 
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Substituting Eq.(6) into (5), yields 

. h dp mt 



dt 



A 



'Hint: p, 



int 



where 



Uintit) = U(t ,t)Hi(t)U(t,t ). 
The magnetic moment operator and the magnetization become 

M = M + Mi 

and 



M 



where 



and 



V 
Mo 



Tr{p int (M ,int + M IM t)}, 
SH 



5H 



ex1 



Mi = -A 



i 



5H 



(8) 

(9) 

(10) 

(11) 
(12) 

(13) 



cxt 



while Mo,int an d Mi^ n t are related with Mo and Mi, respectively, in the same 
way as Hint is related with Hi. The operators Mq, a = 1,2,3, fulfill the SU(2) 
magnetization algebra 



M a o,M{ 



ih 1 e abc M c o, 



(14) 



where the summation over repeated indices is assumed. 

The operators Mo,int, Mi^ n t are generally used to calculate the magnetic suscep- 
tibility [9]. Let us show now how these operators determine the time evolution of 
magnetization. The evolution in time of Mo,int is given by the equation 



dM . 



int 



dt 



h 



U(t ,t) Ho, Mo U(t,t ) 



= lM ,int x H e// . 



(15) 



It describes the magnetization precessional motion with respect to H e jj. 

The equation for Mi^ n t describes more complex magnetization dynamics governed 
by the interaction Hamiltonian Hi- However, this dynamics includes the precessional 
motion as well. Introducing 



D; 



Ho, Mi -iMixHeff 



(16) 
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to represent deviations from the purely precessional motion, we bring the equation 
for M.i t int into the following form 



dMl, int Ca tt 

— = lM IMt x H e// 

+U(t ,t)(^- + B I )u(t,t ). (17) 

Taking the time-derivative of Eq.(ll) and using Eqs.(8),(15) and (17), we finally 
obtain 

dM - 

— = -| 7 |MxH e// + D, (18) 

where 

B = \(Um,Hi] ) + (^ + D 7 ). (19) 



dt 



Therefore, Eq.(18) is the magnetization equation for the system specified by (1). This 
equation is general since it is derived without specifying Hi. The D-term contains all 
effects that the interactions, Hi, can have on the magnetization precession, so that 
Eq.(18) is complete. 

^ — * 

The contribution of Hi to the D-term in the magnetization equation can be di- 
vided into two parts. One is proportional to {[Ai,Hi]_) and is related to the change 
in the density matrix when the interactions of Hi are switched on. The second part 
{^f - + D/) originates from the change in the magnetization itself. Which part of D 
is dominating depends on the nature of the interactions. For the interactions related 
to the relaxation processes, the D-term represents a general form of magnetization 
damping. 



3 Example: spin-spin interactions 

The spin-spin interactions among the spins in the system introduce many body effects, 
which can be treated perturbatively in the weak coupling regime. In this case the D- 
term can be expanded in powers of A. To demonstrate this, we consider the spin-spin 
interactions of a specific type. The interaction between spins is usually an exchange 
interaction of the form 

-2J£^. = -— Ml (20) 

the coupling constant J being called the exchange integral. We generalize the ansatz 
given by Eq.(20) by assuming that the exchange integral depends on the magnetiza- 
tion and introduce the spin-spin interactions as follows 

\Hi = Y.J ab { M )StS^ (21) 
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where J ab = XM a M b . Since Tti does not depend explicitly on the external field, its 
contribution to the magnetic moment operator vanishes, hAi = 0. 

The non-vanishing commutator [A4 ,7ij]_ in Eq.(19) is the only contribution of 
the spin-spin interaction to the magnetization equation, resulting in 



where 
and 



D = -M x ft, 
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Q a = (M a ,M b )M b , 



M a ,M b = M a M b + M b M a 



4b \Aa 



(22) 

(23) 
(24) 



The correlation function G ab = ( M.q, Mq\ ) is the sum of spin correlation functions, 



(25) 



i j^i 



excluding the self-interaction of spins. For the standard ansatz given in Eq.(20), 
D = and the magnetization equation does not change. 

If the spin-spin interactions are turned on at t — t , so that p(to) — Po(to), then, 
integrating both sides of Eq.(8), we find 



P\(t) = Po(to) + ^J t dr 



n x (r),p x (r) 



(26) 



Substituting Eq.(26) into the definition of G ab , yields the equation 

G ab (t) = G% b (t ) 

+ 1 f drJ cd (T) (e ace G ebd {T) + e bce G aed (r)) , 

7 J to v ' 



where 



Ct 



Ml Ml X 



(27) 



(2? 



which relates G ab to the third order correlation function, i.e. the correlation function 
of the product of three magnetic moment operators, 



G abc = ( M a ,M b 



,M C r 



>■ 



(29) 



The correlation function G abc , in turn, is related to the fourth order correlation func- 
tion and etc., and we have therefore an infinite number of coupled equations for spin 
correlation functions. For any practical calculation this infinite hierarchy has to be 
truncated. That then defines the approximation scheme which may be considered 
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on the basis of the physical requirements for the system. The approximation scheme 
will depend on the physical properties such as density and on the strength of the 
interactions. 

If the Hamiltonian Hi is a small perturbation to the original Tio, we can solve 
Eq.(26) perturbatively. In the lowest, zeroth order in A, we replace p\(t) by po(to), 
so that G ab Gg 6 (t )- We choose the initial value for G ab as 

EE(^)o = ^ (so) 

i j^i 

with I xx = py = 0, I zz = I and I ab = for a ^ b. We also define again the z- 
direction as the direction of the effective magnetic field that is chosen uniform and 
static. Then the D-term becomes, in component form, 

D x = 2\ 1 IM y M Z) (31) 
D y = -2X 1 IM X M Z , (32) 
D z = 0. (33) 

producing two effects in the magnetization equation: the magnetization is now pre- 
cessing with respect to (H z — 2X1 M z ), its z-component remaining constant in time, 
(d/dt)M z = 0, and the frequency of the precession is 

cJ = ^o(l-A— — ). (34) 

Therefore, in the lowest order of perturbations, when the D is linear in A, it shifts the 
direction and the frequency of the precessional motion without introducing damping 
effects. To find a role of the higher powers of A in D and to determine how they affect 
the magnetization equation, a truncation of the chain of spin correlation equations 
is needed. This would require a consistent perturbation approach to the hierarchy of 
the coupled equations for the correlation functions. 



4 Conclusion 

We have derived a general form of magnetization equation for a system of spins 
precessing in an effective magnetic field without specifying the internal interactions. 
It can be applied in the study of magnetization dynamics of any type, including non- 
equilibrium and nonlinear effects, provided the interaction of individual spins with 
each other and with other degrees of freedom of the system is specified. 

— * 

The D-term in the magnetization equation has been obtained without using any 
approximation scheme. It is exact, accumulates all effects of the internal interactions 
on the magnetization precessional motion and can be a starting point for practical 
calculations. For the spin-spin interactions, it is determined by the spin correlation 
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functions, which fulfil an infinite chain of equations. A further analysis of the D-term 
requires an approximation scheme to truncate the chain in a consistent approach to 
higher order calculations. 

In our talk, we have considered a specific type of the spin-spin interactions, which 
do not contribute to the algebra of magnetic moment operators. However, if spin- 
spin interactions depend explicitly on the external field, the form of the algebra can 
change. In this case, the total magnetic moment operator becomes nonlinear in Aio, 
and this results in the magnetization algebra with an infinite chain of commutation 
relations. The chain has to be truncated in a way consistent with the truncation of 
the chain of equations for the spin correlation functions in the same approximation 
scheme. 

Numerical computations along the lines developed in [10] can provide a further 
insight into the problem. This talk is based on the work [11] where an extended list 
of references can be found. 
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